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ABSTRACT 

We present a novel description of how energetic electrons may be ejected from the pulsar 
interior into the atmosphere, based on the collective electrostatic oscillations of interior elec- 
trons confined to move parallel to the magnetic field. The size of the interior magnetic field 
influences the interior plasma frequency, via the associated matter density compression. The 
plasma oscillations occur close to the regions of maximum magnetic field curvature, that is, 
close to the magnetic poles where the majority of magnetic flux emerges. Given that these os- 
cillations have a density-dependent maximum amplitude before wave-breaking occurs, such 
waves can eject energetic electrons using only the self-field of the electron population in the 
interior. Moreover, photons emitted by electrons in the bulk of the oscillation can escape along 
the field lines by virtue of the lower opacity there (and the fact that they are emitted predomi- 
nantly in this direction), leading to features in the spectra of pulsars. 
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1 INTRODUCTION AND CONTEXT 

The source of the plasma in pulsar magnetospheres is a subject 
that has long been at the heart of pulsar electrodynamics, and is 
a problem still open to question (for example, see review articles 
by [Miche l ( 1982, 20041). It is assumed that there must be an initial 
electron flux from the stellar surface, and that these particles are 
then implicated in the production of the electron-positron plasma 
which populates the magnetosphere ( [Ruderman| 1 97 1 [ |Sokolov^ 
|Temov|1986[|Harding & Lai|2006^ . The stellar surface in question 
is the transition between the pulsar interior and its atmosphere; the 
formation of an environmental pair plasma then must be linked to 
the energetics of electrons expelled from the outer crust into the 
atmosphere. 

Therefore we are interested primarily in the physical mech- 
anism that extracts electrons from the interior, and projects them 
into the atmosphere immediately above the pulsar surface, where 
the electron-positron plasma is formed. This breaks down the fun- 
damental pulsar problem into two parts: 

(i) From where do the energetic electrons come? 

(ii) Given these energetic electrons, what underlying physical 
processes cause the creation of the atmospheric pair-plasma? 

This article attempts to contribute a new perspective on (i), 
the source and energetics of ejected electrons. The central idea 
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is that electrons are ejected from the magnetic flux tubes as a re- 
sult of the self-field induced in plasma oscillations inside the tubes 
themselves, by virtue of the magnetic curvature. Electrons are as- 
sumed to be severely constrained to follow magnetic field lines in 
the interior; these lines are strongly curved near the poles, where 
the magnetic field emerges normal to the pulsar surface (this is be- 
cause the magnetic field to some extent must run parallel to the 
curved surface in the interior, but then emerge in a small spot that 
constitutes the polar region; all emergent field lines are sharply re- 
fracted as they leave the iron crustal interior - a shell of iron sur- 
rounding a denser neutron core - and emerge into the atmosphere). 
The curvature of the field along which electrons are constrained to 
move causes density variations in the electron population, as parti- 
cles are forced to lose parallel momentum. Such density variations 
cause electrostatic compressions and rarefactions of the electron 
'gas' at the plasma frequency, with the self-electric field providing 
the restoring force (that is, the mutual repulsion of the negatively 
charged electrons). 

These electrostatic waves also have a characteristic radiation 
signature at the plasma frequency, since the electrons are strongly 
accelerated when participating in this bulk motion. Given that re- 
gions of high magnetic flux have a relatively low opacity for pho- 
ton transport, this means that photons radiated in this way are more 
likely to move parallel to the field and emerge at the pole. Given 
that the plasma density in the pulsar atmosphere is substantially 
less than that prevailing in the interior, the atmosphere is transpar- 
ent to these photons, and we present tentative evidence that such 
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signatures are already apparent in the x-ray spectra of selected ob- 
jects. 

Tiie plasma frequency plays a key role here in the electrody- 
namics, and we show how the plasma frequency must be dependent 
on the magnetic field, since the latter causes the matter deformation 
and compression that enables the perfectly conducting, anisotropic 
transport of electrons along flux tubes to exist. 

The sections that follow address key issues in sequence: the 
magnetic compression, and its influence over the plasma frequency; 
the dynamics of 1-D electrostatic waves and oscillations; the mod- 
elling of the Fermi energy, and the associated maximum energy 
gain in an electron wave; and a selection of possible candidate pul- 
sars showing non-thermal features that correlate with the prediction 
of radiation at the internal plasma frequency. Concluding remarks 
finish the paper, and an explanatory Appendix gives more detail on 
the relativistic calculation of the Fermi energy. 



2 DENSITY COMPRESSION AND THE PLASMA 
FREQUENCY 

The process of extracting electrons from the outer crust, requires 
the description of the structure of matter there, and this has to be 
supported with the study of the structure of atoms in strong mag- 
netic fields. 

There are various exotic descriptions of atoms in extreme 
magnetic fields, showing how the basic lattice structure of the con- 
ducting metal is severely distorted, with implications for the con- 
duction electrons. Such models are controversial, and in some re- 
spects, difficult to reconcile (for example, the literature ^Neuhauser] 
|et al.|1986||1987[|Lieb et al.|1992[ > is undecided about the complex 
nature of the bonding between iron atoms under such conditions). 

However, the more general concept is widely accepted, 
namely that the magnetic compression of the atoms leads to an ef- 
fective iron atomic radius R given by jLai|200lt 



7?! 



,1/5,-2/5 



ao 



(1) 



where Z is the atomic number (26 for Fe), b = B/Bq is 
the ratio of the pulsar magnetic field to the critical field Bo ~ 



2.35 X lO^'T and ao = 5.29 x 10"''m is the 



Bohr radius. For the typical pulsar field of 10 T, 6 ~ 426 and 
7i « 9xlO"^^m.TheclassicalradiusofFeis_RFe = 1.4xl0~^"m 
jSlater|1964^ , so that the compression leads to an increase in matter 
density of {Rpe/Rf ~ 3.75 x 10'\ Of course, the free electron 
density is also increased as a direct consequence of this compres- 
sion. Taking the free electron density of Fe under standard terres- 
trial conditions as n^" = 1.7 x lO^^m"^ jAshcroft & Mermin| 
|1976| l, we can arrive at an electron number density for the pulsar 
from atomic compression by applying the same scaling: 



X 3.75 X 10^ 



: 6.4 X lO^V"^ 



(2) 



which is a lower limit, since it is likely that there are more electrons 
able to access the conduction band under compression than in the 
classical limit. 

However, there is another way of deriving the empirical elec- 
tron number density in the pulsar crust. Iron has a monatomic body- 
centred cubic (bcc) crystal structure under standard terres trial con- 
ditions, with inter-atomic spacing apc = 2.87 X 10"^°m (jAshcroftl 
|& Memun|1976|. As suming that the iron on the surface of the pul- 
sar is also bcc ( Ruderman | 1 97 1 | l , this means that the compression 
factor here can be calculated in terms of the effective atom spacing 
in the pulsar environment compared with the terrestrial one. 



The presence of extraordinarily large magnetic fields distorts 
the iron crystal from being isotropic to being severely anisotropic 
in transport terms: the iron is very highly conducting in the direc- 
tion parallel to the magnetic field, but perpendicular transport is 
severely inhibited. The overall picture is of a set of perfectly con- 
ducting 'tubes' aligned locally with the ambient magnetic field di- 
rection along which electrons are able to move relatively freely; this 
description is qualitatively consistent with the 'distorted atoms in 
a strong field' model. This is essentially the structure proposed by 
|Ruderman| ( |1971^ ; |Canuto & VentiiraH1977^ show that the conduc- 
tivity parallel to the magnetic field is on average 20 times greater 
than that of the field-free case; conversely, the transverse conduc- 
tivity has a strong dependence on the magnetic field and the Fermi 
energy, and is typically orders of magnitude smaller than the longi- 
tudinal case. 

With this simple picture we can capture the essence of the 
electron motion; parallel momentum is unconstrained, but per- 
pendicular momentum is quantised. In Ruderman's simple one- 
dimensional flux tube (Ruderman 1971 ), the radius of the flux tube 
is set equal to the mean orbital radius associated with the Landau 
ground state: 



P = 



h 

eB 



1/2 



2.6 X 10^ 



(3) 



where we have taken B — 10*T. 

In so doing, the simple I-D flux tube model assumes domi- 
nant parallel motion along the magnetic field, and assumes that the 
transport anisotropy is sufficient to render all non-trivial Landau 
levels as unimportant, effectively confining the electron to a single 
flux tube. 

The scale-length p is then the effective inter-atom spacing 
in the compressed bcc structure, leading to a compression of 
{oF^ / pf' ~ 1.4 X 10*", and an electron number density n^L based 
on lattice compression of 



: X 1.4 X 10** 



2.4 X 10^^ m" 



(4) 



Since we now have possible electron number densities ranging over 
more than 2 orders of magnitude, we shall take the geometric mean 
as the characteristic pulsar interior electron number density n^: 



Fc ,6/5 
XAO 
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2.61 



= n^XL^'^ « 2.4 X 10^" m"^ 
XL = aUeBo/nf'^ ^ 159 
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Jie = (n^n^) =n^[XAXL) o 
« 3.5 X 10^" X fci-^'^ m-^ 



1.23 X lO^'^m"^ 



(5) 

(6) 
(7) 

(8) 



where we have assumed a magnetic field of _B = 10*T, that is, 
b ^ 426. 

The geometric mean compression yields a mass density of ap- 
proximately 10''kgm~^ for ^'^Fe, which is consistent with the con- 
ventional assumptions of density on the surface ( [Shapiro & TeukoP] 
|sky|1983||Harding & Lai|2006l l. 

The electron number density is a critical parameter for plas- 
mas, since it defines the basic collective timescale, namely the 
plasma frequency uip. Given that the interior is an excellent exam- 
ple of fixed positive ions and mobile electrons, we can define the 
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plasma frequency here in terms of the electrons only: 

1/2 



1.05 X lO^'fe"-^"^ 



(9) 
(10) 



where we have used the geometric mean number density in the nu- 
merical evaluation; me and eo are the electron rest mass and electric 
constant, respectively. 

Notice that the plasma frequency in the pulsar interior depends 
on the magnetic field strength; it is not an independent parameter. 
This is an important result that, although appearing to be a straight- 
forward deduction from the mass density dependence well-known 
in the literature, is still worth stressing because of the radiative con- 
sequences. 

Having used a blend of the published equations of state for the 
bound material in the lattice to infer the free electron density, it is 
important to recognise that such free electrons are not required to 
obey the same equation of state as the bound matter: after all, the 
electrons in a cold plasma have no equation of state! Since the cold 
electron plasma is an excellent model for free electron behaviour 
in our pulsar interior, we shall therefore not require any further de- 
tailed analysis of the interior equation of state. 



3 1-D ELECTRON DYNAMICS 
3.1 Calculating the Fermi energy 

The simple assumption of such a strongly anisotropic conducting 
flux tube has major implications. The electron dynamics in the 
interior can be treated essentially as one-dimensional, since the 
electron momentum parallel to the magnetic field greatly exceeds 
that perpendicular to it; hence in modelling the electron distribu- 
tions, and deriving the Fermi energy as a substitute for the sur- 
face work function, a one-dimensional treatment will be a good ap- 
proximation. Moreover, since cross-field transport is inhibited, the 
role of the Landau levels in the interior is diminished: in the non- 
relativistic formulation of the electron motion in a uniform strong 
magnetic field B jSokolov & Temov|1986| l, the electron energy e 
is given by 



+ pf|/(2me) + (n + \)hLUc 



(11) 



where u)c = eB /rae is the cyclotron frequency, is the electron 
rest mass and n = 0, 1, 2, ... is the principal quantum number for 
the quantised electron orbit in the plane perpendicular to the mag- 
netic field. The radius R of orbit of the perpendicular motion of 
the electron can be expressed (in the same non-relativistic limit) as 
ISokolov & Temovl ( [T986| l 



R ■- 



eB 



1/2 



(12) 



The mean orbital radius associated with the Landau ground 
state is given by 

/ \ 1/2 

R{n = Q)={ — \ Ri2.6xlO"^V (13) 

in other words, p = R[n — 0), in Eq. ([sjl, as already assumed. 
Note that the energy increment Se between Landau levels is 



of 10 T, Landau levels are separated by ~ 11.5keV, implying that 
there is negligible population of the higher Landau levels if ther- 
mal excitation is the only mechanism, given that the typical surface 
temperature of a pulsar is < lO^K (in energetic terms, < lOOeV) 
( Haberl 2007, Kargaltsev & Pavlov|2007, Bogdanov et al. 200^. 
This reinforces the merit in assuming that the electron momentum 
is largely parallel to the magnetic field, and lends credence to the 
argument that a one-dimensional statistical treatment captures the 
essential physics. 

Finally, although calculations involving the Landau levels are 
quantum in nature, they are not relativistic (though the generalisa- 
tion is possible). For the moment, we will defer the detailed discus- 
sion about the need for a relativistic treatment to the Appendix. 

In general, the distribution function for electrons in the pres- 
ence of a magnetic field B is really the Fermi-Dirac distribution, 

.fpD, 



fpD = {exp[(e ± figB - p)/{kBT)] + 1}' 



(15) 



where e is the energy, /i is the chemical potential, /i_g is the mag- 
netic moment of the electron , and T is the temperature. In the limit 
T ^ 0, fj, ^ ep, the Fermi energy, and since hbB ~ 6keV 3> 
kgT for B ~ lO^T, we can assume that the electrons are spin- 
aligned in the lowest energy configuration, and so this term can be 
neglected. Hence we have 



lim fpD 



^irn{exp[(6-M)/(fcBT)] + l}- 



1 for e < /I 
for e > fi 



(16) 



For temperatures such that fcsT ^ /i, the distribution is therefore 
basically a step function, with all energy levels equally occupied up 
to fi, and none of the higher ones occupied. 

IRuderman] ( | 1 97 1 [ l calculated the Fermi energy for the simple 
1-D case, motivated by the restricted motion of the electrons im- 
posed by the enormous magnetic field strengths in the pulsar inte- 
rior. The reason for calculating the Fermi energy is that ep is an 
excellent guide to the work function of the surface, that is, the po- 
tential barrier which must be surmounted before interior particles 
can escape to the exterior. 

Assuming that the mean electron energy is far below the Fermi 
temperature (so that the step-function nature of the Fermi distribu- 
tion can be assumed), for electrons in the population. 



N : 



9{'^)fFD{e)de : 



(e)de 



(17) 



where g[e) is the density of states for the electron gas, and /i = 
is the Fermi energy in the limit of T ^ Tp- For the one- 
dimensional electron gas, g{k)dk — [L/(27r)]dfc, where L is the 
characteristic scale-length of the problem, and we have suppressed 
the normal degeneracy factor 2, given the spin-alignment assump- 
tion. In this case, the integration yields 



(18) 



where we should interpret N/L as the line-density of electrons. 
Assuming a uniform density approximation within the pulsar inte- 
(14) rior, the line density of electrons confined to a flux- tube of radius p 
is simply the volume electron number density times the cross- 
where B is given in Tesla. Hence for a typical pulsar magnetic field sectional area associated with the flux-tube: N/ L = yrp^ . When 



& = —B ^ IQ-^B eV 



© 0000 RAS, MNRAS 000,[T]{8] 



4 D. A. Diver et al. 



substituted, this gives 

1,2 2 2 'A 



2me 



nl 

8mee2 52 
1.03 X 10"' 



Magnetic field lines 



£2 



2.28 X 10"^*^ 6" '' J 



(19) 



(20) 



in SI units. Eq. \i9) is essentially Ruderman's formula ( [RudermMil 
[T97T] >. 

Assuming a typical magnetic field of lO^T, and taking the 
electron number density from ([8l, yields 



^ 97keV 



(21) 



Notice that the Fermi energy increases with magnetic field 
strength: this is because the electron number density increases 
faster than the field strength, no matter which scaling model is 
chosen. This is an important, if slightly counter-intuitive, result, 
since the influence of the magnetic field in expressions for ep such 
as Eq. l[T9j suggest that an increasing magnetic field strength will 
lower the potential barrier at the poles. However, this neglects the 
indirect influence of B on the electron number density: increasing 
field strengths must lead to increased matter compression, and this 
must be taken into account in the Fermi energy calculation. 

Note that since the surface temperature of the pulsar is be- 
low lO^'K, equivalent in energy terms to < lOOeV, then \2l\ is an 
acceptable approximation, since in this case it is self-evident that 
Tp « Tp. 

The calculations presented here assume a 1-D model for elec- 
tron motion, and therefore preclude any contribution from Landau 
levels. The Appendix details the nature of the calculation where 
Landau levels must be included, and where the energies are intrin- 
sically relativistic. 



3.2 Electron motion in inhomogeneous magnetic fields 

The simple 1-D description of the electron transport requires care- 
ful treatment when the magnetic field is not homogeneous. Mag- 
netic field emerges from the interior predominantly at the poles; 
this concentration of the field lines at specific points on the sur- 
face means that there must be significant curvature inside (and out- 
side) the star, as the field emerges from inside the iron spherical 
shell. Electrons constrained in the 1-D model to follow these field 
lines are nevertheless inertial particles, and therefore cannot instan- 
taneously change direction. However, they are also restricted in 
their perpendicular motion, since they have to satisfy the Landau 
quantization. In the simple 1-D model, the electrons negotiate this 
magnetic field curvature in a non-quantal way, since scattering into 
Landau levels above the ground state is not permitted; therefore 
the electrons must change direction to follow the field direction, 
without any significant locally transverse excursions; momentum 
conservation must be provided by the crystal lattice, and the over- 
all energy conservation must lead to electrons losing energy to the 
crystal as they negotiate the curve. 

The overall effect of the magnetic inhomogeneity must be to 
introduce a local 'bunching' of electron density along the magnetic 
field direction in the region of greatest local curvature, that is, near 
the magnetic poles. Note that complex magnetic geometries can be 
produced by thermal evolution of the neutron star itself ( Page 2004 , 
|Pons|2009l l, where strong toroidal field curvature can result from 
the strong magnetic feedback on thermodynamic transport proper- 
ties. If strongly curved magnetic fields produced in such process 




Figure 1. Schematic showing a simplified sketch of part of the magnetic 
field of the pulsar, interior and exterior (note that only a few magnetic field 
lines are included, in order to keep the diagram simple). A small region 
near the pole is magnified to show the location of the electron oscillations 
induced by the curvature; the shaded regions are intended to illustrate den- 
sity compressions. Note that strong magnetic field curvature can also be 
caused by magneto-convective processes not represented in this diagram. 



lead to significant perpendicular field eruption at the pulsar sur- 
face, then these field lines will also act as electron ejection sources. 
A schematic of the situation is shown in Figure 1 . Such local den- 
sity fluctuations will drive an electrostatic wave along the flux tube, 
accelerating (and decelerating) non-local electrons (ahead and be- 
hind the curved region) that populate this tube. Now either this is 
sufficient in itself to eject electrons at the end of the tube, where the 
surface exists, or the electrostatic wave itself evolves non-linearly 
and eventually breaks, leading to the ejection of a few relatively en- 
ergetic particles ahead of the wave. Either way, it seems that such 
a description has the requisite element of a parallel acceleration 
mechanism that is self-consistent, and not dependent on any frame- 
transformed field component. 

It is worth distinguishing between the electron temperature 
and the electron energy; mono-energetic beams of electrons are for- 
mally cold (that is, possess zero temperature), yet can stream in 
the beam direction with significant energy. In the scenario outlined 
above, it is assumed that the mean energy of the electrons, includ- 
ing any directed streaming, is less than the Fermi energy, since oth- 
erwise a significant fraction of the interior electrons would simply 
leak out of the flux tube into the pulsar atmosphere at the surface. 
Whilst this might solve the electron production problem, there re- 
mains the issue of identifying the physical process that causes such 
streaming. For the purposes of this article, therefore, we will as- 
sume that the electrons are mostly trapped in the iron crust, just 
as the free electrons are trapped in a metal in a natural terrestrial 
context. As in the latter case, there needs to be a mechanism that 
provides the necessary energy for electrons to overcome the po- 
tential barrier at the interface between the pulsar surface and the 
vacuum (or atmosphere), before electrons can be liberated at the 
surface. Thermionic emission seems implausible, since the surface 
temperature is less than one tenth of the Fermi temperature; hence 
our concentration on the (self) electric field produced by charge 
concentrations of electrons in density waves. 

On a more formal basis, electrons with energy less than the 
escape energy (that is, the Fermi energy) that are directed towards 
the surface along the curving magnetic field lines must be reflected. 



© 0000 RAS, MNRAS 000,[T]j8] 



On the surface extraction of electrons in a pulsar 5 



since they cannot escape. Therefore there must be two populations 
of electrons in motion along a flux tube; those that are moving to- 
wards the surface, and those that are moving towards the interior, 
having been reflected at some point. This is similar to the magnetic 
bottle effect ( |Boyd & Sanderson|2003 ) where charged particles are 
reflected at magnetic field concentrations; note that the analogy is 
not perfect, though, since in the magnetic bottle, particles are re- 
flected simply because energy is transferred from parallel motion 
along the field to perpendicular motion in the form of larmor or- 
bits. Such counter-streaming electrons will lead to density instabil- 
ities if the perturbations on each stream happen to evolve in phase 
such that they reinforce each other. The full dispersion relation for 
electrons moving parallel and anti-parallel to the magnetic field di- 
rection with speed u is given by |Boyd & Sanderson] p003) 



This result can be extended to the relativistic case: 



2 2 
_|_ ^^2 



(22) 



in which — n-ie^ / [eora^) is the square of the plasma frequency 
for electrons moving parallel to the magnetic field with speed u, 
where ni is the number density of such electrons; is similarly 
defined for electrons moving anti-parallel to the field. We have as- 
sumed symmetry in the modelling, for simplicity; clearly we expect 
m « 712, which is equivalent to saying that the electron loss-rate is 
small; hence we also expect u)i ~ ujp, i = 1,2. The frequency and 
wavenumber of the common perturbation on the electron streams 
are given by u) and k respectively. We must solve Eq. l[22j in or- 
der to arrive at the criteria for instability. Manipulation of the al- 
gebra reveals that the dispersion relation is biquadratic in u), and 
that there are imaginary components of frequency (or wavenum- 
ber) when u < ypluj^jk. Put another way, any density perturba- 
tions with a wavelength A satisfying 



1-KU 



<\<R, 



(23) 



must be unstable, where we have taken the upper limit on the wave- 
length to be the maximum scalelength in the star itself, that is, the 
pulsar radius ii, . Since the lower limit is at most 27rc/(\/2i^p) ~ 
5 X 10^ ^"m ~ 200/5 ^ _R,, then it is always possible for instabili- 
ties to affect the counterstreaming flows, in the form of longitudinal 
plasma oscillations at the plasma frequency. 

There are further assumptions in this simple treatment, namely 
that the plasma density is uniform, and the electrons are monoen- 
ergetic. These assumptions are fine in the context of an idealised 
cold plasma treatment; we contend that they serve as an adequate 
illustration in the context of the pulsar flux tube. 

The nature of this streaming instability is to induce electro- 
static oscillations on the dynamic electron population inside the 
flux tube. 

These plasma oscillations produce a self electric field because 
of the compression and rarefaction of electron density; there is no 
associated magnetic perturbation, since the particle and displace- 
ment currents cancel perfectly. Note that this is a pure electron 
plasma with a fixed positive ion background, and so density en- 
hancements here are pure electron overdensities. There is however 
a maximum amplitude of electric field associated with such oscil- 
lations. If the wavenumber is k, then in the non-relativistic limit, 
this maximum electric field iSmax l |Kruer|1990[|Mori & Katsouleas| 
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where 7,^ = (1 + p|/(m^c^))^''^ is the relativistic factor for max- 
imum associated electron momentum in the oscillation. This is 
the limiting field before coherent motion breaks down and wave- 
breaking occurs, at which point part of the electron population be- 
gins free-streaming at speeds greater than the phase velocity of the 
oscillation itself ( [Mori & Katsouleas|19"90l|Kruer|1990| l. 

Consider the case of an electrostatic oscillation in a flux tube, 
very close to the surface near the magnetic pole. If the oscillation 
is at the point of wavebreaking, then the maximum energy 5e that 
can be gained by an electron falling through the self-field of the 
oscillation before free-streaming must be given by 

5e ~ eEmax x distance travelled by electron 

This electron will escape to the surface if it has an energy greater 
than the Fermi energy. Hence for escaping electrons. 



eEmaxT^/k > €f 



which can be expressed as 



k' < 



5.6X lO^V-^^m-^ 



(26) 



(27) 



(24) 



using Eqs. |4]l, (jTO| and \20\ . For the typical pulsar magnetic field, 
b ~ 426, yielding k < 1.7 x lO^^m"^. In other words, plasma 
oscillations with wavelengths (or scale variations) A — 2n/k ex- 
ceeding 4 X 10~^^m can eject electrons from the surface, providing 
such oscillations are at maximum amplitude. This minimum wave- 
length is considerably greater than p ~ 2.6 x 10^^^, the width of 
the flux-tube (and the effective inter-ion spacing), suggesting that 
whilst the breakdown of ultra-short wavelength oscillations can't 
provide sufficient energy for electrons to escape the surface, longer 
wavelength oscillations will readily do so. 

Thus we have shown that instabilities in the electron dynamics 
within the flux tube are sufficient to eject electrons from the pulsar 
interior. 

Given that the opacity of the pulsar to photons of frequency u) 
is reduced by a factor (cj/loc)^ for propagation parallel to the mag- 
netic field direction ( |Lai|200l||Canuto|I975[ l, where luc = eB/m 
is the electron cyclotron frequency, there is the enhanced prospect 
of photons characteristic of the oscillation within the flux tube es- 
caping directly to the pulsar atmosphere; in this context, lj — u)p, 
and {ujp/iUcf ~ 6Ab-° '^^ « 0.12 for magnetic fields of 10**T, 
leading to the enhanced probability that the magnetic polar regions 
are transparent to the emission of keV photons (from the plasma os- 
cillation) produced inside the flux tubes by the bulk plasma dynam- 
ics. Note that for magnetic fields of around 4 x lO^T, this opacity 
enhancement disappears, and so direct emission of photons below 
0.4keV by this process wouldn't be expected. 

There has been considerable discussion of the electromag- 
netic transmission properties of the pulsar crust, particularly in the 
strong-field limit relevant to magnetars, in which the free electrons 
in the crust influence the transmission of the blackbody radiation 
from the pulsar surface itself (see for example [van Adelsberg|(2006l 
|2005[ >; [Turona et al.|f2004^ ). Strictly, this article presents a discus- 
sion of the electrostatic (not electromagnetic) behaviour of elec- 
trons oscillating at the interior plasma frequency. These electrons 
radiate in the local oscillating electric fields close to the surface 
such that some of the radiated photons escape from the interior. 
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and are able to propagate in the pulsar environment; a fraction of 
them may survive atmospheric processes to reach the observer. It is 
helpful therefore to clarify that the characteristic non-thermal sig- 
nature of such a process as described in this article is not affected 
by the cold plasma twin-mode electromagnetic absorption proper- 
ties of the crust as discussed in the literature. 

Of course, the picture is further complicated by the fact that 
the self-field of the oscillation can accelerate particles in the per- 
pendicular direction ( [MacLachlan et aT]|2009 V Whilst the 1-D 
model forbids perpendicular transport within the flux tube inside 
the pulsar, this restriction is no longer valid at the transition to the 
atmosphere (since the confining crystal lattice is no longer impos- 
ing dynamical constraints), and it is possible that some of the emit- 
ted electrons will possess a perpendicular velocity component rela- 
tive to the magnetic field direction very close to the surface, leading 
to a spread in exit trajectories for escaping electrons. 



4 SUPPORTING FEATURES IN NEUTRON STAR X-RAY 
SPECTRA 

We offer the following examples from XMM-Newton observations 
in support of our hypothesis that there is non-thermal radiation at 
the internal plasma frequency present in the pulsar emissions. The 
motivation for searching for such a non-thermal signature in X-ray 
spectra is that its presence in an object's spectrum lends plausibility 
to the electron production process actually occuring in that object. 

It is established in the literature that X-ray pulsars commonly 
have a significant soft X-ray excess l |Hickox et al.|2004) l, but not 
all such excess signals come from accretion processes (see for ex- 
ample points 3-5 in the conclusions of [Hickox et al.|[2004| l. We 
offer here some examples in which we suggest that certain X-ray 
specral features can be associated with the novel physical processes 
described in this article. Sifting all available data in support of non- 
thermal emission that could be attributed to plasma frequency radi- 
ation from the poles is a daunting task, not just because such sig- 
nals are magnetic field and geometry dependent, and may in fact lie 
outside the spectral range for XMM-Newton. However, we did find 
several examples, which are described below. The magnetic field 
streng ths in each example (other than in the case oFDe Luca et al.| 
j2004j ) were taken from the loffe pulsar catalogue[^] 

Please note that the predicted plasma frequencies are pulsar 
surface values; terrestrial observation of such features will incur a 
gravitational redshift of up to 25%. 

De Luca and co-workers ( |De Luca et al. 1'2004) showed de- 
tailed observations of the x-ray spectrum of the neutron star IE 
1207.4-5209, principally to demonstrate the presence of cyclotron 
absorption features at 0.7, 1.4, 2.1 and 2.8 keV, from which was 
derived a surface magnetic field strength of B ~ 6 x lO^'T, signifi- 
cantly different from the field of 2.6 x 10**T derived conventionally 
from timing parameters. The interpretation of the higher harmonics 
has been challenged l |Mori et al.|2005) as an instrumental effect. 
However, assuming that the lower magnetic field strength is cor- 
rect, and using the plasma frequency given by Eq. for this field 
strength, we calculate ujp ~ 9.4 x lO^^rad s~^, generating photons 
with energy hvp of approximately 0.62keV. This is very close to an 
excess above blackbody radiation in the phase integrated emission 
in the band 0.1-2.5 keV, just before the absorption feature at 0.7keV 
(see Figure 6 in |De Luca et al.|2004j . Caution must be exercised 



http://www.ioffe.nj/astrol/psr-catalog/Catalog.php 



here, since there is also a sharp drop in the CCD efficiences in both 
the pn and MOSl instruments at 0.5eV ^Kirsch„2004j ; however, 
since the observations show a rise above the blackbody spectrum 
before 5eV, followed by a sharp drop, it is reasonable to assume 
that the excess is genuinely observed. 

PSR J0538-I-28I7 ( |McGowan eFaLKIOOSl [Zavlin & Pavlov] 



2004 1 has a surface magnetic field strength of 7.33 x lO'^T, giv- 
ing hvp « 3.3keV. The energy spectrum of this pulsar shows a 
marked excess over the blackbody emission at 3 keV (see Figure 
4 in M cGowan et al.|2003 1 consistent with non-thermal emission 
from the pole via the mechanism reported here. There are no cali- 
bration issues at these energies with the PN and MOS I instruments 
on XMM-Newton, both of which recorded the excess signsil. 

The surface magnetic field of Geminga (Zavlin & Pavlovj 

2004 1 is 2 X 10*T, implying that radiation from the sub-surface 
plasma oscillation will be at energies ~ 6.6keV. The fitted EPIC 
count-rate spectrum shows a clear departure from the model black- 
body spectrum in the range 6 — 8 keV, as shown in Figure 6 of 
[Zavlin & Pa vlov ( 2004). 

PSR J1932-I-1059 (B 1929-1-10) shows a smafl feature at 2keV 
jWozna et al.| (j2003j. Figure 4) which is not fitted by either the 
blackbody spectrum or the power law. The pulsar has a surface field 
of — 5.18 X lO^T, suggesting that the plasma frequency yields 
photons of around 2.6 keV, offering a plausible explanation of this 
spectral feature. 

If the identified spectral features are associated with the 
plasma frequency, then this means that there are indeed plasma os- 
cillations present at the poles, and therefore it is also possible that 
such oscillations are contributing to the production of high-energy 
electrons in the atmosphere. Whilst it is arguable that the features 
we are citing here are close to being borderline in terms of statis- 
tical significance, we contend nevertheless that they offer a tanta- 
lising glimpse of what might be the seat of the electron expulsion 
mechanism, and as such, merit attention. 



5 CONCLUDING DISCUSSION 

The general scenario presented in this paper is as follows: plasma 
oscillations are induced in the throat of the tightly-curving mag- 
netic field structure as it exits the pulsar interior at the magnetic 
poles. Such oscillations can sustain maximum amplitude electric 
fields such that at wavebreaking, electrons can be accelerated to 
energies in excess of the Fermi energy and ejected from the inte- 
rior into the immediate pulsar environment. This behaviour could 
be intermittent, since the transit time for electrons to move along 
a magnetic field line from one pole to the other is less than a typi- 
cal pulsar period (less than a ms, in fact), and so there could be an 
influence on the sub-pulse structure. 

In addition to ejecting energetic electrons, the accelerations 
associated with these interior plasma oscillations will give rise to 
radiated photons which escape along the low opacity field lines 
and emerge at the magnetic poles. Since the electron density im- 
mediately above the magnetic pole (produced by the ejection of 
electrons from a lossy mirror) is much less than that in the interior, 
any plasma oscillation radiation that leaks out will propagate away 
- photons of around lO'^'^Hz (around 2 keV) could be directly emit- 
ted from the polar region by interior oscillations - another source 
of moderate energy radiation. We have presented evidence of X-ray 
spectra of pulsars which show behaviour of this kind. 

Given that the self-field in the interior of the flux tube is 
the source of the acceleration, this field will not be uniformly 
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aligned with the magnetic field, leading to the possibility of highly- 
energetic imperfectly aligned electrons very close to the pulsar sur- 
face. Such electrons could gyro-rotate (subject to Landau levels) 
and radiate, or perhaps even collide with iron nuclei; either way, 
energetic photons are possible directly from the surface. 

Of course, the pulsar still loses electrons this way; charge bal- 
ance must be maintained ultimately, and this is probably achieved 
by a combination of sucking in electrons from the electron cloud 
above the poles (each half-cycle the self-field of the oscillation will 
reverse) and surface diffusion of electrons back into the interior 
from the exterior. Note that any positive potential arising from elec- 
tron depletion of the pulsar will increase the work function (as dis- 
tinct from the Fermi energy), making the potential barrier higher for 
electrons to escape the surface. It is conceivable that such charging 
considerations could be sufficient to suppress electron ejection at 
the surface, even to the point of influencing the pair-production in 
the atmosphere. 

This general picture has several favourable aspects: not only 
does it provide a plausible source of parallel acceleration, leading to 
electron production at the magnetic poles, but also there is the real 
possibility of producing azimuthal structure in the electron produc- 
tion at the surface, arising from the mutual interaction of waves 
in adjacent flux tubes. Moreover, this model is consistent with the 
ideas associated with the distortion of atoms under strong magnetic 
fields giving rise to sheaths of free electrons aligned with those 
fields; it also yields accelerating parallel electric fields, albeit from 
a different physical process. 

It is also worth noting that dropping the magnetic field 
strength raises the Fermi energy at the surface (since the strong field 
is suppressing the work function); this in turn reduces the spectrum 
of electrostatic oscillations that can eject free electrons, reducing 
the efficiency of atmospheric electron production and the associ- 
ated energetic beaming processes. This may give an insight into 
transient pulsar behaviour, for those objects with weaker surface 
fields 
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7 APPENDIX 

7.1 Calculation of Fermi Energy in a Landau Levels 
environment 

The general case should take into account the effect of Landau lev- 
els on the Fermi energy. In fact the Fermi energy is dependent on 
the linear density of the electrons along the magnetic induction field 
lines, and the regime becomes relativistic above certain values. Fol- 
lowing this idea the energy e — jmc^ is given by 



7= < 



spm up 



pI eB 
IH + 2{n + l)h—- spin down 



(28) 



Then we see that the same value of 7 has two possible values (n 
and n -\- 1) with the exception n = 0. Then the total number of 
electrons is given by 



N 



211% 



E 



dp\\ 



exp 



kT 



(29) 



+ 1 



because the transverse moment is filled in Landau levels. 

The value of the Fermi energy is ej? = limr^o M and there- 
fore 1 16 1 implies, using py — zu^^mc^, that there are upper limits 
for n and ccz. Putting £ = e~m(? = eimc^and/i — fiimc^ , k — 
fiLo/mc^ with 



1 + ■Tom -f nn - 



(30) 



The upper limits of ci7|| are reached when yL\ — ei and there- 
fore the limits of integration are 



whereas in the limit 



N h 
2L mc 



L mc ', 



Ml 



(39) 



N h L mc 

2Lmc'^ '^2]VX 
Eq. is similar to Ruderman's ( |Ruderman| 1 97 1 ^ expres- 
sion for Fermi energy when we neglect the spin effects and the zero 
order energy. This means it can only be applied to subrelativistic 
regimes. However the general solution Eq. l |37[ ) also contains the 
relativistic case for high density electrons or large magnetic field 
amplitudes. 

The overall solution shows that the value of the Fermi energy 
depends on the value of the magnetic field, and if 



2L mc 



< 



heB 



(40) 



it may even be the dominant term determining its value. However 
for the particular case of exclusively longitudinal motion, there 
should not be any dependence on the magnetic field; that it enters 
the Fermi energy even in this case is puzzling. 

Moreover the combination of a quantum statistical description 
in the perpendicular plane and a classical one along the magnetic 
field leads to the situation in which the statistics may preclude oc- 
cupation of any Landau levels other than the zero one, as given in 
Eq. l[35j. These contradictions between the description in Landau 
Levels and the underlying statistics have to be resolved in future. 



+ 2/i — 2nK, 



spm up 



/ii + 2/i — 2{n + 1)k spin down 



0, mt{ 



+ 2^1 , 



After some algebra we find that 



iV h 
L mc 



int[{M?+2Mi)/«;] 

E 1 




flf + 2/il — UK 



+ \hi'i + 2/il - {n + l)K 



(31) 
(32) 

(33) 
(34) 



an equation in /ii whose solution has to be found numerically, pro- 
vided 



^Il + 2fii 



> 1. 



(35) 



However if only n = is present, this equation turns into 
N h _ 
L mc 
with general solution 



(if + 2/ii + + 2Ati — 1^ 



Ml 



. N h L mc 

^ ' 2Lmc^'^2Nli' 



1/2 



i N h L mc , , , ^ , 
If ( h K < 1, then Eq. d37b reduces to 

2L mc 2N h 



N h L mc 

2 \2Lmc ^'^2iVX 



(36) 



(37) 



(38) 
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